Y-system for form factors at strong coupling in AdS§ 
and with multi-operator insertions in AdS% 



Zhiquan Gao° and Gang Yang c 



O 

■ a State Key Laboratory of Theoretical Physics, 
$_i . Institute of Theoretical Physics, Chinese Academy of Sciences, 

j^i Beijing 100190, China^ 

b Kavli Institute for Theoretical Physics China, 
Beijing 100190, China 

c II. Institut fur Theoretische Physik, Universitat Hamburg 
Luruper Chaussee 149, D-22761 Hamburg, German^ 

> 
00 

; Abstract 

^ ■ We study form factors in N =4 SYM at strong coupling in general kinematics and with 

multi-operator insertions by using gauge/string duality and integrability techniques. This 
generalizes the A0IS3 results of Maldacena and Zhiboedov in two non-trivial aspects. The 
first generalization to AdS^ space was motivated by its potential connection to strong 
coupling Higgs-to-three-gluons amplitudes in QCD which was observed recently at week 
coupling. The second generalization to multi-operator insertions was motivated as a step 
towards applying on-shell techniques to compute correlation functions at strong coupling. 
In this picture, each operator is associated to a monodromy condition on the cusp so- 
lutions. We construct Y-systems for both cases. The Y-functions are related to the 
spacetime (cross) ratios. Their WKB approximations based on a rational function P(z) 
are also studied. We focus on the short operators, while the prescription is hopefully also 
applicable for more general operators. 
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1 Introduction 

One of the most challenging problems of modern theoretical physics is to understand 
the dynamics of strong coupling QCD analytically. While this is still very difficult, lots 
of studies have been focused on simpler models such as theories with supersymmetry. 
The general philosophy is that a good knowledge of these theories may finally help us to 
understand the real QCD. A particular interesting theory that has drawn much attention 
is the J\f—4 super Yang-Mills theory. There have been a number of evidences that J\f—4 
SYM results are important building blocks of QCD quantities, see for example [JJI2]- By 
the gauge/string duality, it becomes also possible to study the M—A SYM in the strong 
coupling regiem where it is dual to a perturbative or semi-classical string theory in an 
AdS background 011115]. 

An impressive achievement is that we now have a good control of computing anomalous 
dimensions in M = 4 SYM to any reasonable order in principle, see for example [6l[71[8l[9j 
[TO] , where the intergrability of the theory plays a fundamental role [TT | [T2 | [T3] (for a review 
on many aspects of integrability see [H]). It is expected that the similar achievement may 
also be made for other more complicated observables, such as scattering amplitudes and 
correlation functions. Indeed surprising dualities and integrable structures have been 
found for amplitudes and null Wilson loops [T^[16l[T7l[T8l[T9l|20l|2TJ , and also correlation 
functions in a special light-like limit [22|[23] . 

One remarkable development is the computation of scattering amplitudes in Af = 4 
SYM at strong coupling [T5] . It was shown that the problem can be dual to a string 
minimal surface problem in AdS. The solving of this non-trivial geometrical problem 
was developed based on the integrability of the classical worldsheet theory p4 l l25 | l26]R 
Hopefully, these classical results will be useful to solve the full quantum problem such as 
in the study of operator dimensions [3"2"]l3"31. 

In this paper, we will focus on a more general class of observables, the so-called form 
factors. They are observables involving both on-shell particles and off-shell operators, 

1 See also [27] for a treatment of tricky 4if-gluon cases, [28,29 for the study of Regge limit, and [30|l3T] 
for the connection to CFT in the regular-polygon limit. 
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therefore are in some sense hybrids of amplitudes and correlations functions 



(Out-states | Y\ ^A x i) I In-states) . 

i 

We will consider form factors in pure momentum space 

F(qi,--- ,qi;Pir-- ,Pn) = II f d A x l e i ^(0(x 1 )---0(x l )\p 1 ---p n ), (1.1) 

fc=i •* 

where pf = and qu is arbitrary. 

Most studies so far have been focused on the form factors with one operator inserted. 
Form factors in string theory in AdS were first studied in [34J. Based on the recent 
developments of strong coupling amplitudes, a T-dual picture of form factors was proposed 
in [55], and the problem was solved in the AdS 3 case by using integrability techniques 
in [36]. At weak coupling, form factors in Af = 4 SYM were first studied in [37J, and 

have received attention only recently, see for example [38j[39j|30lllTJ|l2] ' ^ ne P ar ^ cmar ^y 
surprising observation in [43] is that the remainder function of a two-loop three-point 
form factor in = 4 SYM matches exactly with the maximally transcendental part of 
the two- loop Higgs-to-3-gluon amplitudes in QCD [44"] Fl. 

As this correspondence looks very intriguing, one may think that this is an accidental 
coincidence. However, this two-loop coincidence is already rather non-trivial, which may 
be appreciated by a simple look at the very different perturbative structures of Feynman 
diagrams in M = 4 SYM and QCD. It may be therefore reasonable to expect that there 
could be some hidden relations which will explain this coincidence and might play further 
roles for other situations, at least for the three-point case due to its particularly simple 
kinematics^. If the two-loop coincidence is going to be true for higher loops, one may 
expect the strong coupling form factors would carry a non-trivial piece of information 
of strong coupling QCD. Considering that there are very few tools to study strong cou- 
pling QCD amplitudes, this possibility provides us enough motivation to study the strong 
coupling form factors seriously. 

The computation of form factors at strong coupling in [36] was restricted to two 
dimensional kinematics. In this case the non-trivial quantities start at four-point. In order 
to study the three-point form factor, one needs to consider more general kinematics. In 
this paper we will consider the form factors in full R 1 ' 3 kinematics, corresponding to string 
in AdS§. As is usually happened, the generalization from AdS% to AdS$ is a nontrivial step. 
Although the underline picture is similar to the AdS^ case, the monodromy structure in 
AdSc, is more complicated. In particular, the truncation conditions involve small solution 
contractions which are not T-functions. This complexity also makes it much more difficult 



2 The relation between form factors and Higgs-gluons amplitudes may be understood by noting that 
the operator in form factor [43] is equivalent to the Higgs-gluon effective vertex obtained by integrating 
out a quark loop. 

3 It would therefore be interesting to study three-loop case. Hopefully the progress can be made in 
Af = 4 side, as in [33] (also with the techniques developed in [IS]), while the computation in QCD seems 
much more challenging. 
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to construct the Y-system, which is a main new challenge of the AdS§ problem. We will 
describe the general construction, and the Y-system for the three-point form factor will 
be explicitly given. 

Another interesting generalization in this paper is to compute form factors with multi- 
operator inserted. The main motivation is to study correlations functions at strong cou- 
pling with the help of on-shell techniques. Similar idea has been used at weak coupling 
in [IB]. Although the observables we consider contain on-shell structures, they involve 
multiple operators, and in principle should contain all kinds of information of correlation 
functions. In particular, one should be able to extract the OPE coefficients from form 
factors containing two or more operators. 

The basic idea we propose may be illustrated by the following flow chart@ 



Cusps 



Small solutions 



Hirota system 



Y system 



Operators 



Monodromy matrices 



The main picture is that for each operator one can define a corresponding monodromy 
matrix, which will give a linear relation for the small solutions. These small solutions are 
related to the cusps and are the same building blocks for calculating amplitudes, therefore 
the known method of computing amplitudes can be applied to these more general class of 
observables. It is in this sense that we can compute off-shell observables by using on-shell 
techniques. 

We derive explicitly the Y-system for form factors with multi-operator insertions in 
AdS 3 , while in principle it should be possible to generalize to the AdS 5 case. The construc- 
tion proposed in this paper is expected to be in principle applicable to arbitrary operators, 
although the study will be focused on light operators!!, for which the monodromy can be 
given explicitly. 

This paper is organized as follows. In Section |2} we review the main physical pictures 
and the general strategy of strong coupling computation via AdS/CFT and integrability. 
We then review form factors in the AdS$ case in section [3j Form factors in general 
AdS$ kinematics are developed in section [U and the three-point case is discussed more 
explicitly in section [51 The generalization to multi-operator insertions is given in section [61 
In section [71 the P(z) function and WKB approximation are studied. Section [8] contains 
a summary and some discussions. There are three appendices. Appendix [A] is a collection 
of the definition of T and Y functions and their corresponding equations. A review of 
(momentum) twistor variables is given in appendix [B] Appendix [C] is a brief discussion of 
the monodromy in a different basis. 



One may note that this is different from the logic used in computing anomalous dimensions via 
Y-system. Here it is important to obtain the Y-system, where the Y-functions are interpreted as the 
spacetime cross ratios, and for which the boundary condition can be conveniently introduced. 

5 These include short BPS operators such as the stress tensor supermultiplets which are also mostly 
studied for form factors at weak coupling, and also light non-portected operators with dimensions oc A 1 / 4 . 
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2 Classical string and integrable system 



Due to the nature of the problem which involves a few different stories and intermediate 
steps, in this section we give a brief review of the whole picture. The discussion here is 
not supposed to be self-contained, but we hope to cover the key physical pictures and 
central ideas. We suggest interested readers to the original papers (in particular [2UES]) 
for more details. 



2.1 Form factor as a classical string solution 

As a first step to set up the problem, we explain how to map the observables at strong 
coupling in the M = 4 SYM to a classical string problem in an AdS background. Our 
focus is on amplitudes and form factors developed in [T5J[35]- 

We first consider the picture for gluon states. Recall the AdS space in Poincare 
coordinate 

ds2 = dy^dy + d_l (2l) 



z 



Gluon states in = 4 SYM are dual to open strings on the IR D3 branes (as an IR 
regulator) at the horizon (i.e. z — > oo) [15]. One important property of the open strings 
on IR branes is that they carry very large proper momenta. Because the high energy 
scattering is dominated by a saddle point approximation [47], the computation of open 
string amplitude becomes a classical string problem. 

Form factors also contain operators, which are dual to closed string states in the bulk 
with boundary condition at z — > [HE]. Therefore form factors correspond to open and 
closed strings scattering coming from the horizon and the boundary respectively, as shown 
on the left-hand side of Figure llR 

To simplify the problem, one important trick is to apply a formal T-duality along y^ 
directions [48j[T5]0. The T-dual space is still an AdS space 

2 _ dx^dx^ + dr 2 

where r = 1/z. The boundary and horizon reverse their roles in the T-dual space. 
The momenta of strings become the "windings" of strings. For amplitudes the problem 
becomes a type of Wilson loop problem J51]|52] , but with a null polygonal boundary. For 
form factors, the boundary becomes a periodic null Wilson line [35], where the period is 
determined by the momentum of the closed string q. The minimal surface also extends 
to the horizon, as illustrated on the right-hand side of Figure [XFI . 

Therefore, the form factor problem becomes to find the area of the minimal surface 
(over one period) with boundary conditions at both the boundary and the horizon of the 



6 It is assumed that the scattering is still dominated by the classical saddle point. 

7 This is in the sense of using Buscher's formalism defined at action level [49 , in which it is also 
straightforward to generalize to fermioinic directions [50] . 

8 It is also obvious that a mixing of Wilson loop and operators such as studied in [S3] are very different 
from the form factors we considered here. 
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boundary horizon boundary horizon 

Figure 1: The picture of T-duality for form factor, q is the momentum of the operator 
which corresponds to a closed string state in the bulk, q = J2iPi due to momentum 
conservation. After T-duality, the picture becomes a minimal surface ending on a periodic 
null Wilson line at the boundary and extending to horizon. The period is given by q. 



T-dual AdS space. The general structure of the strong coupling results is 

Observable = e~"^~ Area x (string a' corrections) , a' ~ —= . (2-3) 

v A 

The string corrections in principle may be computed by considering string fluctuations 
where the classical solution is taking as a background, along the line of [54j|j. 



2.2 String in AdS as a classical integrable system 

Because of the non-trivial boundary conditions, it is very hard to solve the string equa- 
tions. The idea, proposed in [21] (see also for example [SHUSHIED]), is that rather than 
solving the string equations directly, one can apply the Pohlmeyer's reduction ]§T\ to 
reformulate the string equations to a Hitchin like system and then use the techniques of 
integrability. Here we review this main strategy. 

Since Pohlmeyer reduction is a well-understood procedure, we only point out that 
after the reduction, the string equations of motion and Virasoro constraints take a form 
of flat equation 

dA- z -dA z + [A z ,A z ] = 0. (2.4) 
If one decomposes A into two parts A = A + <3>, the equations form a Hitchin like system 

D z $ g = , D g $ x = , [D z , D z ] + [$ z , $,] = , (2.5) 

9 It seems no such computation has been done for any solutions corresponding to amplitudes, even for 
the simplest four-point case, where both the classical solution [15] and the result (given by ABDK/BDS 
ansatz [55,. 55]) are known. The pure spinor formalism [57] might be useful for such computations. 
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where D z := d z + [A z , ]. For general AdS 5 case it is SU(A) systemH while for AdS% it is 
reduced to SU(2). The flat connection is not arbitrary but satisfies a Z± automorphism 

A = -CA T C-\ <t> z = -iC<5> T z C~\ <&- z = iCQZC- 1 , (2.6) 

where C is a constant matrix whose explicit form is not important here. This constraint 
plays an important role in the construction as we will see later. One can solve the linear 
equation 

(d + A)ip = 0, (2.7) 

where the solution if) is related to the target space coordinates and therefore to the string 
solutions. 

A natural logic would be to first find the solution for A which solves the Hitchin 
equations, and then solve the linear problem to find the solution ip which gives the classical 
string solution. However, the strategy used here is different. Roughly speaking, we will 
use the properties of the linear solution and the flat connection to construct the area 
directly without knowing the explicit solution. 

The key idea is to use integrability. The integrability can be understood by the fact 
that one can lift the flat connection to a family of connections 

A -> A{() = (A z + ^Adz + {A g + C^z)dz, (2.8) 

while the Hitchin equations are still satisfied. The new parameter ( is called spectral 
parameter. We also use another variable 9 where ( = e %e . If one solves the linear problem 
with .4.(0, one obtains a one-parameter family of solutions ip(C), and the original physical 
solution can be obtained by taking ( — 1. 

With this extra parameter it seems one is dealing with a more general problem. How- 
ever, new powerful techniques are available based on this new parameter!"]. The main 
result is that a set of functional equations, so-called Y-system, can be constructed. The 
non-trivial part of the area can be extracted from the solution of Y-system, which turns 
out to be the free energy in a thermodynamic Bethe ansatz (TBA) form [011115]. For 
amplitudes in AdS 5 it is like 



A iree = zZ^J ^ cosh e log [(i + y 1)S )(i + r 3)S )(i + ^ s )^] • (2.9) 

s 

In this form, the area is a function of mass parameters which are implicitly related to 
the physical cross ratios. An important observation later in [66J is that the area can be 
also written as a critical value of Yang- Yang functional, and in the new form the area is 
expressed directly function of cross ratios. 



10 Here we have changed to the spinor representation of 50(2,4). This change of representation is 
equivalent to the using of momentum twistor variables at weak coupling [62) . 

n We would like to point out the idea of introducing new parameters has played many other important 
roles in theoretical physics, such as the 57-deformation in localization techniques and the orbifold gener- 
alization in ABJM theory, see for example the talk by John Schwarz [63 . It would be very interesting to 
study their possible connection to integrability. 
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2.3 Boundary condition and function P(z) 

In this section we explain an important aspect of the story: how to introduce the boundary 
conditions. This will involve a very important holomorphic function P(z). We also discuss 
the special feature of form factors in which an operator is inserted. 

One particular equation of the Hitchin system is the generalized sinh-Gorden equation 

dda -e a - e~ a \P(z)\ = , (2.10) 

where P and a are given as 

P = d 2 X ■ d 2 X , a = \og(dX ■ dX) . (2.11) 

P(z) is a holomorphic function. By making a field redefinition and introducing a new 
coordinate w by a worldsheet conformal transformation as 

a(z,z) = a{z,z) + ^\ogP(z)P(z) , dw = (P(z)) 1/4 dz , (2.12) 

one can simplify the generalized sinh-Gordon equation as 

8 W 8* a - (e & + e~ & ) = 0, (2.13) 

which is a simple sinh-Gordon equation. One should note that the change of worldsheet 
coordinate is only well-defined locally. 

An important fact is that the four-cusp solution (first found in [15]) is simply the 
solution P = 1 and a = a = of the generalized sinh-Gordon equation [60]. This is an 
important reason of doing the above transformation!^!. Asymptotically, the solution near 
each cusp should be same as the four-cusp solution. This implies that near boundary 
where z — > oo, we should have a — > 0. It also implies that each cover of w-plane contains 
four cusps. Therefore, P(z) should be a polynomial, and the degree of the polynomial 
would depend on the number of cusps. The corresponding picture is shown in Figure El 
The coefficients in the polynomial would encode the shape of the polygon. 

12 This makes it also simpler to introduce cut-off and compute regularized area in w-plane |24[|25j . 
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Figure 3: The picture for form factor. An insertion of a operator introduces a singularity 
in the z-plane, and corresponds to a multi-cover of z-plane. 



A new feature of form factors is that there are also operators. As observed in [36], an 
insertion of an operator will introduce a pole term in P(z). This requires a study of the 
boundary condition near the horizon, which will be discussed and generalized to AdS$ in 
section [7J Due to the insertion of operators, the z-plane is no longer smooth. One can 
however smooth the z-plane with the sacrifice of introducing a multi-branch-cover of z- 
plane, as illustrated in Figure [30 This picture is consistent with the periodic null Wilson 
line picture in target space. We will use this picture later to introduce the monodromy 
for small solutions. 

2.4 Small solutions 

Now we introduce the most important building blocks, the so-called small solutions. Con- 
sider again the linear problem 

(d + A{(,z))iP{(,z) = 0. (2.14) 

Because of the special null-cusp boundary conditions, the solution has different asymptotic 
behaviors near different cusps. Small solutions are the solutions which decay fastest while 
approaching the boundary. The small solutions are unique up to a normalization. At first 
sight, it may be confusing why it is the small solution rather than the big solution that 
is important. This can be understood by considering the AdSs case (similar picture also 
applies for AdS^). While approaching an edge i, the solution can be approximated as 

V> ~ cf g ^ + cr al V (2.15) 

It is not the big solution, but the coefficient of the big solution that contains the boundary 
information. This coefficient can be extracted by contracting the full solution with the 
small solution c/ s ~ (ip,Si). In this way, all non-trivial boundary information can be 
obtained in terms of the contraction of small solutions. 

The coefficient c^ ls is related to target-space variable, the momentum twistor Aj in the 
general AdS$ case. A" carries spacetime indices a. On the other hand, the small solution 

13 This is in some sense similar to what happened in w-plane in the 4if-cusp case [24][27] . 
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sf is a solution of worldsheet theory and carry internal-space indices a. This change of 
variables from target space coordinates to worldsheet solutions plays a very important 
role in the strong coupling story. Hopefully, there would be some corresponding picture 
at weak coupling. 

We now consider the relation between small solutions and the spectral parameter. One 
important fact is that: the change of the phase of the spectral parameter corresponds to 
the rotation of small solutions (i.e. cusps) 0- The Z4 automorphism mentioned before 
plays a very important role. For example, the Z2 automorphism in Ad S3 gives the relation 
s i+i(C) s i( el7T ()- The contraction of small solutions can be defined as T-function and 

Y-functions, see Appendix [A] Using Z 4 property and some other identities, it is possible 
to construct a finite set of difference equations between this functions. 

While a Y-system is basically a set of algebraic equations given by a set of determinant 
identities, one needs to provide further information such as the asymptotic behavior of the 
corresponding functions, so that the obtained solution is corresponding to the observables 
being studied. This important information can be obtained by a WKB approximation, 
in the limits of the spectral parameter: \(\ — > or 00. In such limits, the contraction of 
small solutions is dominated by an integral along WKB lines that connect different small 
solutions, for example in Ad^o 

( Si , Sj )\^ ~ exp(^^) ~ * ^ • (2-16) 

The integrand ^/pdz/C, is obtained as the dominant term of flat connection ( 12. 8 j) in the 
limit Q — > 0. The WKB lines can be obtained as the parametric curves z(t) which solve 
the equation Im(i s/p/Q = 0. Therefore they are determined by the function P(z) which 
is related to the boundary conditions. These will be discussed further in section [7J 

One can see that the problem is set up as a Riemann-Hilbert problem (see for ex- 
ample [67]): finding the exact functions from their discontinuities (provided by Y-system 
equations) and asymptotic behaviors (WKB approximations). In this paper we will con- 
struct the Y-systems for form factors in AdS§ and with multi-operator insertions. We 
also study the WKB approximations. The explicit solving of the obtained systems will 
be left to a future study. 



2.5 Conventions 

The basic definitions of T/ Y-functions and how to obtain the basic Hirota and Y-system 
equations are summarized in Appendix |A] For the reader who is not familiar with the 
definitions it may be necessary to have a look at the appendix before reading the following 
sections. Below we mention a few important relations and conventions. 

14 This implies an intriguing correspondence between the worldsheet z-plane and spectral £-plane. 

15 The reason that there is a path connecting different small solutions can be understood that when 
computing the contraction one needs to bring the small solutions to a same point in the z-plane. In the 
limit C 0, the solution Si(z) is determined by the integrand J* z , where the subindex i should be 
understood as the point where the cusp lives. 
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We will often assume the normalization conditions [26] 

AdS 3 case : (si,s i+1 ) = 1, (2-17) 
AdS 5 case : (s;, s i+ i, s i+2 , s i+3 ) = 1, (2.18) 

unless indicated otherwise. The Z 4 automorphism impose the following relations [26] 

AdS 3 case : s<+i(C) = ^3S;(e i7r C), (2.19) 

AdS B case: s m (C) = C^e^C) , s i+1 (C) = C^e^C) , (2.20) 

where 

Si := Si_i A 5^ A s i+ i . (2.21) 

<73 and C are some constant matrices whose explicit forms are not important in this paper. 
There are two different conventions used for AdS 3 and AdS 5 : 

AdS 3 case : /± := /(e ±l l C), / W := /(e^C) , (2.22) 
AdS 5 case : /± := /(e^i C) , f [k] ■= f(e^C) • (2.23) 

Since the number of cusps is always even in the AdS 3 case, for convenience we define 

h := n/2, (2.24) 

where n is the number of cusps. 



3 Review of form factors in A<iS3 

In this section we review the Y-system for form factors in AdS 3 [36] . The construction in 
this case is relatively simple, but the basic physical in the later generalizations picture is 
similar. 



3.1 A look at amplitudes 

We first look at the case of scattering amplitudes. For amplitudes, the correspondence 
minimal surface is smooth. The small solutions are single valued on the z plane 

Sj (e i2 *zX) = Sj(z,() ■ (3.1) 

By definition Sj + n is the small solution in the same sector as Sj but after going around 
the complex z plane once. Because they are the solutions in the same sector and the flat 
connection is single valued, they should be proportional to each other 

s J+ n(e l27T z,C) oc Sj (z,0- (3.2) 

This may be also understood from the periodic condition. Note that an arbitrary propor- 
tional constant is allowed. 
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Figure 4: Small solution and monodromy. 



To do the contraction of small solutions, one needs to bring two small solutions to the 
same worldsheet point. Using (13 .11) and (13.21) . one gets that 

s i+h (z,C) oc Si(z,C), (3.3) 

which implies (sj, Sj+n) = 0, or equivalently T^_! = 0. This provides a natural truncation 
for Hirota equations and the Y-system. The corresponding Y-system is given in terms of 
h — 3 Y-function: Y m , m — 1, . . . , h — 3 [25] . 

3.2 Operator as a monodromy 

Now we consider form factors. Since there is an operator inserted, the worldsheet is not 
smooth but contains a singular point. The small solutions are therefore no longer single 
valued on z plane. In other words, they change their value after going around the complex 
z plane, or more exactly, going around the singular point where the operator is inserted, 
as shown in Figure HI 

This effect can be characterized by introducing a monodromy matrix. One can firstly 
choose two linearly independent small solutions as a basis. To be explicit, one can choose 
{sq, Si}. The monodromy is defined as a 2 by 2 matrix fi(C) satisfying 

^(ze 2 -,C) = fi(C)(^)(^C). (3-4) 

Using the Z2 automorphism relation (12.191) . this also fixes the monodromy relations for 
other small solutions. By taking the wedge of the small solutions, one can obtain 

det[fi(C)] = I- (3.5) 

The exact property of Q is determined by the corresponding operator, which can be taken 
as an input of the system. 

By definition, as discussed for amplitudes above, Sj + n is the small solution in the same 
sector as Sj but after going around the complex z plane once. For the same reason, their 
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proportional relation does not change: Sj + n(e l2 ' K z) oc Sj(z). We introduce the proportional 
constant B(Q so that 

s fl (zX) = B(()s (ze- 2m ,(). (3.6) 

By Z 2 relation (|2.19p . this also determines the proportional constants for other small 
solutions, in particular 

SA+ i(z,C) = BW(() Sl (ze- 2 *\0. (3.7) 

Taking the wedge of small solutions and using the normalization condition one can get a 
constraint on B: 

B- 1 = B [2] . (3.8) 
At the same point of worldsheet, one can obtain 

= bco-^co (*)(*. o, (3.9) 

where the proportional constants are written into a diagonal matrix 

B(() = ( 5[2 (C) B ° {0 ^ , det[B(C)] = I- (3.10) 

Unlike amplitudes, s, and s i+ h are not proportional to each other and (sj, s i+ n) ^ 0, so 
the truncation becomes much more non-trivial for form factors. 



3.3 Truncation and Y-system for form factors 

We consider now how to use the above monodromy relation to truncation of Hirota equa- 
tions, and then how to write them into a Y-system, following [36] . 

Firstly, a useful relation is from the trace of Q, which is conformal invariant. From 
(13. 9p one can obtain 

Tr[fi(C)] = B(C)(s , sn+i)(0 - B-\C)(si, s«)(C) , (3.11) 

where we have used the fact that for a 2 x 2 unitary matrix (det[£T| = 1), Tr[Q] = Tr[0 -1 ]. 
This can be further written as 

Tr[n(c)] := T>[n(c)] = B(C)r ft (0-B- l (c)r^ 2 (c), (3.12) 

where C = e~ i( - A+1)7T/2 . 

One can see that the trace relation provides a truncation for the Hirota equations, 
since one can solve for in terms of T^_ 2 and Tr[f2]. As mentioned before, Tr[0] can be 
taken as an input of the system. 

While Hirota equations are not gauge invariant, it is necessary to write the system in 
a conformal invariant way, i.e. to find a Y-system. This can be done by defining a new 
Y-function 

Y(C) := B-\<;)Tm(Q. (3.13) 
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Figure 5: Lattice picture for the Hirota and Y-system of AdS^ form factors. Black dots 
are for functions T a ^ m . Dots with black circle correspond to Y-functions Y a m , while dot 
with red circle is for Y . Red dot is for which can be solved via the trace relation. For 
the amplitudes, the Y-functions are truncated to the left hand side of the red dashed line. 



Then using ( 13. 12[) . T^(Q can be solved as 

2k(C) = B-\C) [Tr[0](C) +F(C)] , (3.14) 

and furthermore for Y^-i 

Y h -i = T h . 2 T h = Tr[fi]F + F 2 . (3.15) 

A equation for Y is also straightforward to obtain which is simply Y + Y = 1 + Yft_2- 
In this way, one obtains a set of equations in terms of n — 1 Y-functions: 

Y+Y- = (1 + ^+0(1 + ^-1)^ (3.16) 
^-2^-2 = (l + Fn- 3 )(l + Tr[fi]F + F 2 ), (3.17) 
Y + Y~ = l + y ft _ 2 , (3.18) 

where s — 1, . . . , n — 3. A lattice structure of the T and Y functions is shown in Figure El 
Comparing to amplitudes, the system involves two new functions Yn-2 and Y. This 
matches with the counting of the degrees of freedom. The new operator is associated to the 
2x2 monodromy matrix. Because det[fi] = 1, there are three independent components, 
which correspond to three functions: Tr[fi], Y^_2 and Y, while Tr[fi] is taken as an input 
in the Y system^. These n — 1 Y-functions also match the number of degrees of the 
T-dual Wilson line picture, which is 2(n — 1). 

Because of the simple structure of the above Y-system, it can be easily written in 
terms of integrable equations. The free energy part of the area can be extracted from the 



16 This is not necessary to be true. In some special cases the trace of the monodromy may also depend 
on the spectral parameter, then one can introduce another Y function while the trace does not appears 
in the Y-system, see section 5.4 of [36j . In the periodic case that we will consider in this paper, the trace 
of monodromy will be fixed to be a pure number. 
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solution and takes a TBA form [36] 1^1. 

-4free = ^~ J d0 cosh9 log{l + Y s ) + 2^- J d9 cosh 6 \og(l + Y) . (3.19) 

For the case here, the construction of Y-system looks straightforward and simple. 
As we will see later, the generalizations to AdS^ and to multi-operator insertions are much 
more involving. However, the underline pictures are similar. We would like to emphasize 
again how the picture of the operator enters the story. 

The basic building blocks are small solutions, which are defined by the null cusps. The 
operator is introduced as a monodromy matrix which imposes some linear relations on the 
small solutions. Therefore, the operator appears only implicitly and does not play a direct 
dynamic role. The method of computing null Wilson loops can then be applied in almost 
the same way to compute form factors. The detailed information of the monodromy is 
taken as an input. In principle one could consider arbitrary operators, depending on the 
choice of corresponding monodromies. 

3.4 Spacetime picture 

Finally, in this subsection we review the spacetime picture [36]. We first clarify the 
difference between the worldsheet monodromy and spacetime monodromy. Then we solve 
the monodromy and also write the Y functions in terms of target-space variables which 
specify the spacetime boundary configuration. 

There are two different definitions of monodromy. One is defined in terms of small 
solutions, as discussed above. It characterizes the non-single valued behavior of small 
solutions while going around the path which surrounds the operator. We call it worldsheet 
monodromy and denote it as Q. The other definition is defined in terms of spacetime 
variables. It will be called spacetime monodromy and denoted by Cl. 

The spacetime monodromy can be taken as a spacetime conformal transformation. It 
can be given by mapping {X n , X n+ i, X n+2 } to {X , Xi, Xo}. Using spinor decomposition 
Xi = Xf Af , it is enough to consider left-hand spinor Af l 18 l. Explicitly, the monodromy 
can be defined 

Af + „ oc flAf = (H l ±H 2 ). (3-20) 

where i = 0,1,2. Because Xi is the embedding coordinate, the mapping is only projec- 
tively, and an arbitrary proportional constant is allowed for each relation. One also has 
det(O) = 1, because the conformal group is SL(2,R). 

The trace truncation equation A3. 11 j) in terms of spacetime variables can be written 

as 

Tr[6](A£,Af> = (A L ,fiAf) + (fiA L ,Af). (3.21) 

17 Here for simplicity m s ,m are choose to be real. 

18 Note that in the AdS$ case, A is used for momentum twistor. 
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Unlike the worldsheet picture where monodromy relates different small solutions, in 
spacetime the monodromy matrix operates on a single twistor variable. The worldsheet 
monodromy fl^ carries the indices of small solution basis, while the spacetime monodromy 
Cl a b carries spacetime indices. For each component of t2y and Q a b, they are in general 
different from each other. But for special combinations such as the trace Tr[fi], the 
worldsheet and spacetime monodromies take the same value [36J. 

One can solve for Cl more explicitly in terms of Porincare coordinate x. It is convenient 
to consider light-cone coordinate x :— Xq ± Xi, which is related to the spinor X L as [21] 

The monodromy relations ( 13.201) written in terms of x + are 

1 fin + Cl 12 xf 

X t+n ^21 + &22%f 



(3.23) 



where i = 0,1,2. Together with the condition det(f2) = 1, they uniquely fix the mon- 
odromy (up to a whole sign). 

To be more explicitly, we focus on the short operators which is T-dual to a null Wilson 
line boundary condition. One has 

4+n = 4 + 1, (3-24) 
for all %. Using H3 . 231) . one can obtain the monodromy 

O = ( I ?V (3.25) 



1 1 

and Tr[f2] = 2. As shown in [36], the trace of the worldsheet monodromy takes the same 
value 

Tr[fi] = Tr[f2] = 2. (3.26) 

The Y function can be also written in terms of spacetime variables. To do this one 
needs first to write the Y function in a form independent of the normalization. One can 
define 3i(z) = si(ze~ l2n ) and therefore {sa, §x) = B(sn, Sn+i)- Recall the definition (13.131) . 
one obtains 

Y(Q = l£ili4( Ce — (n+D/2)_ ( 3.27) 

[Sh, Si) 

This form makes it obvious that the WKB approximation of Y forms a closed path which 
contains the singular point corresponding to the inserted operator [36]. 
Then Y can be written in terms of spacetime coordinates 

nc = = -^Mt = - „. 5^:,^,, = 4^ ■ p.28) 
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One can see that Y is scale invariant, but different from other Y-functions which are usual 
conformal cross ratios. This is related to the fact that form factor is not dual conformal 
invariant, unlike scattering amplitudes. The important point is that at strong coupling in 
the worldsheet picture, the integrability techniques are still available. One can deal with 
Y exactly in the same way as for usual cross-ratio Y functions. 



4 Form factors in AdS§ 

In this section, we study form factors in AdS^. The construction will be parallel to the 
AdS$ case, however as we will see, new problems will appear during the construction. 



4.1 Monodromy 



As in AdS 3 , one can first choose four linearly independent small solutions {s_ 2 , s_i, s , Si} 
as a basis for the general solutions of the linear problem. The worldsheet monodromy is 
characterized by a 4 by 4 matrix 0(C), which is defined by the relational 

/ si \ ( si \ 



So 
S-1 
\S-2/ 



>e 2 ™,C) = 0^(C) 



(4.1) 



s 

S-1 

Taking the wedge of small solutions one has det[0(C)] = 1. 

By definition, one has the same proportional relation that Sj +n (e t2n z, C) oc Sj(z,()- 
We choose the proportional constant B(Q that 



~>n+l 



>,C) 



(4.2) 



Using the Z 4 automorphism relations (I2.20p . proportional constants are fixed for all other 
small solutions, in particular 



S n -l(^,C) : 

s n - 2 (z,C) ■ 
There are also extra constraints from (I2.20p : 

bb® = b®bW 



B(e-^C)s-i(e-^z,0, 
fi- 1 ( e -^ /2 C)s_ 2 (e- 27r ^,C). 



B = 



(4.3) 
(4.4) 
(4.5) 

(4.6) 



One obtains that at the same point of the worldsheet Sj and Sj +n are related to each 
other as 

/Sn+l\ / Si \ 



S n -1 
\Sn-2/ 



{z,Q = ^(0-0(0 



s 

S-1 
\S-2/ 



(4.7) 
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Note that compared to the AdS3 case, here we choose fl 1 in the definition for convenience. 
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where the proportional factors are written into a diagonal matrix 

2T 1 := diag{ J B,( J B[ 2 J)- 1 , J B[- 4 l,( J B[- 2 J)- 1 }, det[B] = 1 . (4.8) 

4.2 Truncations of Hirota equations 

Similar to the AdS 3 case, one can apply trace conditions to truncate the Hirota equations. 
The trace of monodromy is a conformal invariant quantity Their spacetime picture 
will be discussed later. We first consider the sing le trace Tr[Q] = J2i=i n u- Usin S 
one can obtain 

Tr[fi] = B n (s_ 2 , s_i, s , s n +i) + #44 (s n -2, s-x, s , Si) 

+ #33 ( s -2, S n _i, So, Si) + $22 ( s -2, S — 1) S n , Si 

) • (4-9) 

This can be further written as 

Tr[Q] = B u T-fJ - B M T^}\ - B 22 (s_ 2 , s_i, si, s n ) + B 33 (s_ 2 , s , s i; s„_i) , (4.10) 

which provides a truncation for the chain of Hirota equations by expressing Ti jn in terms 
of other T functions. One may worry about the other two terms which are not T functions. 
However, they can be expressed in terms of T functions as will be discussed in the next 
subsection. 

To obtain a truncation relation for T 3 n , it is natural to consider the s» variables and 
the corresponding Q, B, wherein 

Tr[m = Bxx T 3 [n l - 23 44 T{^~1\ - B 22 (s\- 2 , s_i, s x , s„) + i3 33 (s_ 2 , s , si, s n _i) . (4.11) 

Note that fl, B are not new but related to Q, B, see Appendix O 

To obtain a truncation relation for T 2jn , one can consider the double-trace (see also [68] ) 

Tr^[fi] := J2 = \ i^M" ~ Tr[^ 2 ]) . (4.12) 

X<i<j<4 

Using (14. 7p and the definition of T functions, one obtains 
Tr®[Q] = BxxB 22 T^ 1] + B 33 BuT^l 

-Bxx B 33 (s_ 2 j So, S n _i, S n+ i) — B 22 Bxi (s_i, Si, S n _2, s n ) 

+Bxx Bxa (s -i) So, s n _2, s n+ i) + B 22 B 33 (s_ 2 , Si, s n _i, s n ) . (4.13) 

As will be discussed in next subsection, all non-T-function small solution contractions can 
be written in terms of T functions. One may also consider further a relation from the 
triple-trace 

Tr (3) [fi] := Q a % Q kk = ~ (Tt[Q} 3 - 3 Tr[Q 2 } Tr[ft] + 2 Ti[Q 3 }) , (4.14) 

i<j<k 

but it can be shown that it is actually an equivalent truncation relation for Tx >n function. 
20 The relation between Ti >m and T^ tTn simply corresponds to changing Si to Sj, or vice verse |26j . 
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4.3 Recursion relations 

In the above truncation relations, several new small solution contractions are involved. 
In this subsection we show that they can be expressed in terms of T functions by using 
recursion relations. A few new functions will be defined for convenience. 

We define the contractions appearing in single trace conditions as R and S functions 

Rl,m := ( s -2; S-i, Si, Srra+2)' m \ Rs,m '■= (S-2> S_i, Si, S TO+ 2}' 1 (4-15) 

S ltm : = (s_ 2 ,s ,si,s m+2 ) [ " ml , S 3>m := (s_ 2 , s Q , Si, s m+2 ) [ ~ m] • (4.16) 

Using the Pliicker relations reviewed in Appendix|Al one can obtain the following recursion 
relations 

r p+ rp T">+ rp 

td _ d- 1,771+1 1 2,m+l Tj r>— J 3,771+1 j 2,m+l /. 1 -\ 

-Kl,m — ^-l.m-l - ^ ~T , -K3,m — ^3,m-l~^ r — , (4- i 'J 

1,771 1,771 3,771 3,771 

T [2] T [2] T [2] T [2] 

c— 3,m ± 2,m—l q q— l,m 2,to—1 /i q\ 

<->l,m — ^ 1,771-1 yq: 7^+ > a %,m — a Z, m -17p~+ 7^+ ■ l 4 - i5 J 

3,771—1 3,771—1 1,771— 1 1,771—1 

Together with the initial condition 

Ri,o — -^3,0 — ^2,1 , Si,o — 2i,i; = T 3jl , (4.19) 

all R and 5 functions can be expressed in terms of T-functions. 

To consider the contractions appearing in the double-trace condition, we first define 

U\ >m := ( s -2> s_i, s m , s m+ 2)' m \ ^3,m := (s_i > So, s m+ i , s m+ s) ' m 2 ^ , (4.20) 
Vj. iTO := (s_2,So,s m+ i,s m+2 ) [ " m] , V 3tm := (s_i, Si, s m+2 , s m+3 ) [ ~ m ~ 21 , (4.21) 

which satisfy 

rp— rp , m+ rp— rp— rp . T"i — 

rr _ 1,771— 1 2,771+1 I" 1,771+1 2,771 r r 3,771 — 1 2,771+1 "I - 3,771+1 2,777 / , or) N 

^1,771 — ^ , ^3,771 — Tf, , 

J- 1,771 3,771 

rp+ /T-l I /J!— /TT+ T"^ T _l_ T" - 

T/ _ 1,771— 1 2,771+1 "I - 1,771+1 2,771 T r _ 3,771 — 1-^2,771+1 "I - 3,771+1 2,771 /, 

^1,771 — 7f, > V 3,m — rp- • 

3,771 1,771 

The four contractions in (14. 13p are then defined as 

W 1>m (0 ■= (s-2,s ,s m+1 ,s m+3 ) [ - m - 1] , Ws, m (0 := (s_ 1 ,s 1 ,s m+2 ,s m+4 ) [ - m - 3] ,(4.24) 
W 2 ,m(C) := (s_i,s ,s m+ i,s m+4 ) [ - m - 2] , W 2) m(C) := (*-2, Si , s m+2 , s m+3 ) [ -' m - 2] ,(4.25) 
which satisfy 

y- u+ -T7 Tt V7 U? -77 T, + 

rjr "l,77l w 3,77l 1,771 1,771 Tjr 3,m w 1,771 3,771 3,771 /. 
^1,771 = — 1 — , ^3,771 = — ' — , (4.26) 

1 2,m J 2,771 

ttt ^3,771+1 ^3,771 — ^2,771^2,771+2 777 ^1,771+1^3,771 — ^2,771^2,771+2 //O^ 

^2,771 — rT, 5 W 2,771 — rT, • i 4 - 2 ') 

J- 2,771+1 J 2,771+1 

Using (I4.22p and (|4.23p . all W functions can be written in terms of T functions. 

The main lesson in this subsection is that any small solution contraction can be ex- 
pressed in terms of T functions, therefore it is enough to focus on the T functions. 
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4.4 Y- system for AdS§ form factors 



While it is straightforward to introduce the trace relations to truncate the Hirota system, 
it is more challenging to obtain a Y-system. 

We find it mostly convenient to introduce three new Y a functions as follows: 

Yi :=A 1 ^ i F 3 := A 3 ^, Y 2 := A 2 T ^~ 2 , (4.28) 

±2,n-\ J-2,n-l 1 l,n-l-L 3,n-l 

where 

D[-n+l] 

M ■= B™, A 3 := (B^)-\ A 2 := . (4.29) 



It is interesting to notice the relations 

A 2 A 2 Af A 3 A 1 A 3 



(4.30) 



AiA 3 A 2 A 2 

which have appeared for the amplitudes of n — AK cases for three particular combinations 
of Y- functions (see page 27 of [26J). 

The Y functions satisfy the following nice equations 



Y 1 Y 3 1 + Y 3 n _ 2 Y\ Y 3 1 + Yi >n _ 2 

Y 2 1 + >2,n-l ' F 2 1 + l2,n-l ' 

Y 2 Y 2 1 + Y"2,n-2 



Y X Y 3 " (l + Y 1)n _i)(l+*3,n-i) ' 



(4.31) 
(4.32) 



Notice that there are functions Y a n _i on the right-hand side of equations. To have a 
closed system, one needs to solve them in terms of other Y functions. This can be done 
by noticing the following relations!^! 

Y a , n _! = (A" 1 T a , n )F a , a = 1,2, 3, (4.33) 

while A' 1 T a , n can be solved directly using the trace conditions discussed in last subsection. 
Since the trace functions of Q are conformal invariant, the trace equations (and therefore 
Y atn -i) are guaranteed to be able to be written in terms of Y functions: Y a and Y a ,m, m = 
1, . . . , n — 2. This will be shown explicitly later in the three-point case. 
The full Y system for form factors in AdS^ can be summarized as 

^a,m^4— a,m (1 "I" Yt,m+l)(l ~\~ ^4— a,m— l) 



Y a +l,rnYa—l,m (1 ~t~ Yi+l,m)(l ~\~ Y(i - \ ,rn) 

(4.35) 



Y a *4-a 1 + Y a , n -2 



Y a+1 Y a _ L (1 + n+l,n-l)(l + ^-1^-1) ' 

where a = 1, 2, 3, m = 1, . . . , n — 2, and Y a n _i can be expressed in terms of other Y 
functions appearing in the equations. Therefore one obtains a closed finite system in 
terms of 3(n — 1) Y-functions, as shown in Figure |6j 



1 This explains also why that we introduce the Y a functions in the above form. 
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Figure 6: Lattice picture for the T/Y-system of Ad S 5 form factors. Dots with black circle 
correspond to Y- functions Y a>m , while dots with red circle is for Y . Red dots are for 
T a , n which can be solved via the trace relations. For the amplitudes, the Y-functions are 
truncated to the left hand side of the red dashed line. 



Comparing to amplitudes where there are 3(n — 5) Y functions, in form factors there 
are 12 more (for n > 4). This can be understood as follows. Form factor contains 
a 4 x 4 unitary monodromy matrix Q which has 15 independent components. Minus 
three traces, 12 degrees of freedom are left. On the other hand, this does not match 
with the spacetime picture of a periodic Wilson line in AdS$, where it only has 3n — 7 
independent degrees of freedom! 22 !. Note this is different from the AdS^ case, where the 
number of F-functions matches with the number of the degrees of freedom. It implies 
that in AdS§, the monodromy matrix of a short operator is not arbitrary but with extra 
four constraints. In practice this is not a problem, since the WKB approximation of Y- 
functions are determined in the same way by the P(z) polynomial. For general operators, 
it would require further information which is not understood yet. 
A proposal for the free energy 

Following the result in AdS% [36], a natural proposal for the free energy is 



A *» = I ^cosh0io g [(i + y liS )(i + y 3 , fl )(i + n, s )^ 

;i + y 1 )(i + y 3 )(i + y 2 )^" , 

where c is an integer factor which may be fixed by studying some simple limits. 



m 

+c — I dd cosh 9 log 
2n 



(4.36) 



22 This can be obtained by a counting of symmetries as 3(n — 5) + 4 + 4 = 3rt — 7. 3(n — 5) is the 
degrees of freedom for an n-cusp Wilson loop, while a periodic n-cusp Wilson line would break 4 special 
conformal transformation symmetries, and there is also an off-shell momentum q which gives the other 4. 
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4.5 Reduction to AdS± and AdS% 



We consider the reduction following the discussion for amplitudes in [26]. The reduction 
to the AdS± case is simply given by taking 



TUC) = T 3>S (C), MO = Y 3 JC). 



(4.37) 



Therefore for form factors in AdS±, there are only two trace conditions of Tr[Q] and 
Tr^ [Q] to consider. This reduction will be used in the three-point case later. 

We consider further to reduce the system to AdS 3 . Besides the relation T l s = T 3 s , 
the linear problem splits into two decoupled problem denoted by left and right problems. 
In an appropriate gauge 



S 2 k 







S2k+1 







'fc+1 



(4.38) 



where s L and s R are the small solution of the left and right AdS$ problems respectively. 
The left and the right problems are related by a rotation in the spectral parameter 



(4.39) 



The small solution contraction in AdS§ is reduced to 

{s 2i ,s 2k+1 ,s 2j ,s 2l+l ) = - (sf , sf) ■ (4.40) 

One can choose a normalization (sf,sf +1 ) = 1, this corresponds to an unusual normal- 
ization (sj, Si+i, Si +2 , Si + 3) = —1 in AdS$. Most equations in AdS$ become identically 
satisfied except for the nodes of T 2)2 k- For these, they reduce to the Hirota equations in 
AdS 3 . 

The monodromy matrix can be decomposed as 



S-l 

so 

\8-2/ 



\ze 



2 iri 







[21 

where using the relation SR :a = s L a , one can get 



Q L (() 



S-l 

so 
\s- 2 / 



(4.41) 



n 



R 



-2] 



(4.42) 



One can check that the three traces relations in AdS§ exactly reduce to the single relation 
in AdS^. In deriving it one needs use the relation 



-(Tr[n h )+Tr[Q R )), 



(4.43) 



while the minus sign is due to the normalization (sj, s i+1 , s i+2 , s i+3 ) = —1. 

The reduction of Y functions can be summarized as in Figure [3 One interesting new 
feature is that the Y a in AdS$ do not reduce directly to the Y in AdS^. 
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Figure 7: The reduction of AdS§ Y -functions to AdS%. Most Y -functions become trivial 
(i.e. -1, or oo as shown in the figure), while half ofY 2 ^ functions are reduced to AdS% 
Y -functions. Interestingly, the three Y a functions do not reduce to Y functions in AdS%. 
Red dots are related only to T^ functions in AdS%. 



4.6 Spacetime picture 

As discussed in the AdS% case, the spacetime monodromy can be understood as a space- 
time conformal transformation. In AdS§ it is convenient to consider the map for twistor 
variables, and in this representation the conformal group is SU(4). An introduction of 
twistor variables is given in Appendix [Bj 

Monodromy matrix can be defined by mapping four twistors {A n _ 2 , A n _i, A n , A n+ i} to 
{A_2, A_i, Ao, Ai} 

\ l+n ex (4.44) 

where i = —2, —1,0, 1. Note the map is only projectively, and an arbitrary proportional 
constant is allowed for each relation. Since Cl e S77(4), one has det(f2) = 1. Like for 
small solutions, we also define 

A- := e abcd \ b \ c \ d , (4.45) 

and the corresponding monodromy Q is defined as Aj +n oc fl Aj, where i = —2, —1, 0, 1. 
The single trace relation written in terms of spacetime variables corresponds to 

Tr[n](A_ 2 ,A-i,Ao,Ai) = (fiA_ 2 , A_i, A , Ai) + (A_ 2 , A , Ai) (4.46) 

+ (A- 2 ) A_i, fiAo, Ai) + (A_2, A_i, Ao, fiAi) , 

and similar from Tr[f2]. For double trace one has 

Tr( 2 )[6](A_ 2 ,A_i,A 0) Ai} = (fiA_2,fiA_i,Ao J Ai) + (A_2,A_i,nAo J fiAi> (4.47) 

+ (fiA_ 2 , A_i, n\ , Ai) + (A_ 2) flA_i, A , ClXj) 
+ (A_2? f2A_i, fiAo, Ai) + (QX-2, A_i, Ao, fiAi) , 

where 

Tr^[Q] := Yl = \ ( Tr ^ 2 - Tr ^ 2 ]) ' ( 448 ) 

l<i<j"<4 
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One can solve for the monodromy in terms of the Lorentz variables. We will focus 
on the short-operator which corresponds to a periodic boundary condition. The transfor- 
mation between twistor variables and Lorentz variables is reviewed in Appendix [Bj The 
derivation of the monodromy will be given in two different ways. 

Firstly one has the relations ( IB. 61) 



A [l A 4] 
A [l A 2] 



xl 



A [2 A 3] 
A [l A 2] 



(4.49) 



where x ± = x 1 ± x°. Without loss of generality, one can choose the periodic direction to 
be along x 1 -direction 



x 



i+n 



X; 



X i+n 



xf + q. 



Using (I4.44p . one can obtain 



X i+n 



X i+n 



i+1 



^2a ^36 A 



a ^b] 
a A 6] 



A [1 A 4] ~^ A [1 A 2] 
A [l A 2] 

A [2 A 3] ~ % 1 A [l A 2] 



^la ^26 A 

for all z. These relations fix the monodromy uniquely as 



A [l A 2] 



(4.50) 

(4.51) 
(4.52) 



n 



(\ 

o 

% q 



0\ 

1 
1 
-iq 1/ 



(4.53) 



There is another simpler way to find the monodromy. From the definition of twistor 
variable (IB.lOp . one can obtain 



Xi oc (Aj , Hi) , X i+n oc (A, , fii - i q ■ A*) 



or equivalently, 



A/4 



h; Cl Xi 



(4.54) 



(4.55) 



where fl is given by the same matrix as above, when q is along x 1 -directiono A pro- 
portional constant hi is introduced explicitly, which plays the same roles as the B factor 
defined for small solutions. Note the relation A n+ j = 6;Aj. 
One obtains that 

Trffi] = Tr[t] = 4 , Tr (2) [Q] = 6 . (4.56) 

As in [32], one can shown that the corresponding traces of worldsheet monodromy take 
the same value, which can be taken as an input of the Y system. 



23 For general q^, the bottom- left 2x2 block of fl in (|4.53[) is replaced by ie a ^q a a- 
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Next we consider to write Y a function in terms of spacetime variables. To do this, 
one needs first to write them in a form independent of the normalization. To make the 
derivation simpler, one may recall the relation Y a = Y a ^ n _i/ (A~ l T a<n ). Since i^, n _i is a 
normal cross ratios, one can focus on (A~ l T 0>n ), which is similar to the AdS^ case. 

Consider first the Y\ case. One can define s_ 2 (z) = s_ 2 (.2e~* 27r ) = B^ 2 ^s n _ 2 (z), and 
therefore (s_2, s n -i> s n , s n+ i) = 5[ _2 l(s n _ 2 , s n -i> s„, s„ + i). One has the correspondence 

(A^T hn )(C) i!- 2 ' gn - 1,gn ' Sn+1 | (Ce-^). (4.57) 

\ s -2, S n _l, S n , S n+ i/ 

Together with the expression of Ki j2 , this gives 

(s_2, S_i, S n _i, S n ) (s n -2, S n -i, S n , S n +i) 

which looks almost like a cross ratio. It also makes it obvious that the WKB lines of Y 
form a closed contour which contains the singular point corresponding to the operator. 

The normalization-independent form can be written directly in terms of spacetime 
coordinates 

Y /£ _ — (^-2) Ki-ii ^«) ^n+i) (A-i, A n _2, A n _i, A n ) (459) 

(A_2, A_i, A n _i, A n ) (A n _2 ; A n _i, A n , A„+i) 

The other two Y functions can be obtained in the same way. The normalization indepen- 
dent forms are 



T7 / >n _ ( s -3) s -2, S„) (s_ 2 , S_i, S , S n -l) _ i7r H^2 

(S_ 2 , S_i, S n _i, S„) (S_ 3 , S_ 2 , S_i, S ) 

y 2 (C) = | 5 - 2 ' ^ * w+l) ^ (Ce-^) , (4.60) 

(S_2, So, S n )(s_i, S n _i, S n , S n+ i) 

and in terms of spacetime coordinates they are 

(f2A_3, f2A_2, f2A_i, A n )(A_2, A_i, Aq, A n _i) 



n-2 , 



^(C = ^~ 



(A_2, A_i, A n _i, A n ) (A_3, A_2) A_i, Aq) 



t? f/ . in^\ (^A_ 2 , fiA_i, A„, A n+ i)(A_i, A , A n _i, A n ) 

F2(C = e > = ~7a A X A \(\ A W r - ( ^ 

\A_2, A_i, Aq, A n ;\A_i, A„_i, A n , A„ + i/ 

Using the results in Appendix [Bj it is also easy to write them in terms of Lorentz variables. 
The expressions for the three-point case will be given in next section. 



5 Three-point form factor 

In this section we study more explicitly the three-point case. This case is interesting 
because of its potential connection to QCD quantities as reviewed in the introduction. It 



25 



also provides an example which shows explicitly how to write the truncation relations in 
terms of only Y-functions. 

The three trace equations for the n = 3 case are given as 

Tr[Q] = B u Tj 3 ] + B 33 - B 22 R+ x , (5.1) 

Trp = B 11 rjf] + B 33 r 3 7 1 -B 2a i2+ 1 , (5.2) 

Tr^ffi] = B 11 B 22 Tl 2 l + B 22 B 33 T 1 , 2 + B 11 B M Tlj-B u B 33 wl 2 l (5.3) 
where, using the relations in section | 



= ^ + Bsi = ^ I y + r » , (5 .4) 

J 1,1 + 3,1 

W hl = + ——+—z • (5.5) 

J 2,l i 2,lJ 3 ,l J 3,l 

These provide the truncation for the Hirota system by expressing T a ^ 3 through T 0j i and 
T a , 2 - 

To construct the Y-system, one can notice that using the definition of Y-functions, the 
T- functions can be solved in terms of Y-functions (with also A factors defined in ( 14.29ft ) : 

T 11 = -^, T 31 = -^, T 21 = ^ 2 _ , (5.6) 

12,\1 1 r 2,ll 3 r l, 1^3,1^2 

rp ^-2 A 2 A^3 /_ -x 

1 1.2 — j -1 3,2 — j J- 2,2 — • ( 0. f I 

13,1* 2 2 -'2,1-' 1'3 

One can then substitute these expressions for T functions into the trace equations. 
The only thing which may cause trouble is the A factors. They need to be cancelled 
since the expressions should be gauge invariant. Indeed, after a little calculation, one can 
express the trace condition explicitly in terms of only Y-functions: 

Tt^M = 2 ' 2 I I 3 3 2 > x %1 (5 1Q) 

if W F^yW^Y^yW^ 1 j 

T7+ V+ v [3] V [3] V [2] ( Y [2] V [2] ' 

r 3 1 2,1 *3 J 2,l r 2 f y[2] v [2] * 1,1*3,1 

^+^131 r l,l r 3,l 



Ff ] F 2 yjlfyw yMH ^i 1 . 



All A factors are cancelled exactly. This provides a non-trivial consistency check for our 
construction. As claimed before, Y a , 2 can be expressed in terms of other Y-functions and 
the trace functions. The final Y-system is a closed system in terms of six Y-functions: 
three Y a \ and three Y a . 
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5.1 Reduction to AdS± 



A three-cusp periodic Wilson line can be always embedded in an AdS± subspace of AdS$. 
Therefore, one can simplify the system further to AdS^. As mentioned before, in AdS^ 
one has Y\ m = Ys m . The Y-system equations are given as 



(5.11) 
(5.12) 





1 + Y lj2 


y- y+ 
I 2,1 I 2,1 


i + y 2i2 




1 + Y 2A ' 


y2 






1 + Y 1A 


^2 


i + n,i 


Y 2 


1 + Y 2;2 ' 




(i+n, 2 ) 2 



where by the trace condition and also using (I4.56p . 



Y 



2:2 



Y 1 



Y, 



4 - 



Y\ Y 2tl 



y[-4]TFh4] 
2 2,1 1 1 



y- y+ y- 

2 1,1 2 1,1 2 1,1 



y[-2]T7t-2] 
2 1,1 2 2 



+ 



IT^" ^l 



2 2,1 



F^F? Y 2 Y [ - 2] 7f Y\f Yl)1 Y 1>lY ® 



(5.13) 



(5.14) 



T7+ V+ T7 



Y\ Y 21 Y 1 Y 2jl 

2 1,1 2 1 2 2 1,1 



2 1 2 2 



Yo 



Y x Yl 



Y 2 
2 1,1 



Y 2 
2 1,1 



2 2,1 2 2,1 



This is the Y-system which has potential connection to strong coupling leading tran- 



scendental pieced of Higgs-to-3-gluons amplitudes in QCD. Note that one can also use 
fl5.lip - fl5.12p to rewrite them into other forms, in particular to change the phase shift of 
some Y functions. 

The WKB approximation of Y-functions is determined only by the P(z) function 



P(z) 



1 



,2 ' 



(5.15) 



z z* 

which will be discussed in more details in section [7J The degrees of freedom also match: 
the complex number a_i provides two real parameters, while the three-cusp periodic 
Wilson line has also two independent ratios variables. 

The equations (I5.13P and (I5.14p looks a little complicated. In particular, a new fea- 
ture is that some functions have large phase-shift which is beyond the physical strip 
(— 7r/4, 7r/4). This will make it a little more complicated to write them in the form of 
integral equations, as the extra pole contributions need to be carefully considered. We 
will leave this problem to another study. 

Finally, we consider to express the Y-functions in terms of spacetime coordinates. As 
in the weak coupling, it is convenient to consider following variables 



u 



Pl2 



v 



P23 



w 



q 2 



(5.16) 



24 In the sense of first taking a summation of the perturbative leading transcendental results which is 
then evaluated at the strong coupling saddle point. 
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where Pij := Pi + pj. There are only two independent variables since 

I 2 = P12+P23+P3D u + v + w = l. (5.17) 
The Y functions in terms of these variables can be obtained as (see Appendix [B]) 

Fl « = ^ = I ( a 1,A 1' A2 ; A u = 7771 — ~TTT • ^ 

0-2 (A-2, A_i, A 2 , A3) 1/(1 -w) + l 

Y (£ — {) — (Ai, A 2 , A 3 , A4)(A_i, A , A 2 , A3) _ v . 
0-2 6-1 (A_2, A_i, A , A 3 )(A_i, A 2 , A 3 , A 4 ) «w' 

v / > -\ (A-2, A_i, A , A 3 )<A_i, A , Ai, A 2 ) uw 

*J - 77 : — : — t ttt t t — rr — > {O.ZU) 

Ao, A 2 , A 3 ; \A_ 2 , A_i, Ao, Ai/ f 

3 i7r/4\ (A- 2 , A_i, A 2 , A 3 )(A_i, Aq, Ai, A 2 ) w ( 1 



2 ' llC " j " (A_ 1 ,A 1 ,A 2 ,A 3 )(A_ 2 ,A_ 1 ,A ,A 2 ) " l-^ll- W + V' 

One can compare them with the interesting set of variables necessarily appearing at 
weak coupling [13] in constructing functions via the so-called symbol technique [69J: 

1 1 1 1 1 1 uv vw wu\ 
u, v, w, 1 — u, 1 — v, 1 — w, 1 ,1 ,1 , , , > . (5. 22) 

U V W W U V J 

One can see that similar combinations appear in Y functions. This is the same as in the 
six-gluon amplitude case, where the variables in the symbol construction [69] correspond 
to the Y-functions at strong coupling^. The three-point form factor provides a further 
evidence that the "correct" variables for constructing functions from symbols at weak 
coupling, which is hard to know (usually only through guess work), may be read directly 
from Y-functions. 



6 Form factors with multi-operator insertions 

In this section we consider the form factors with multi-operator insertions 

F(qi,--- ,qr,Pi,--- ,Pn) = J] / d i x l e i ^(0(x 1 )---0(x l )\p 1 ---p n ). (6.1) 

k=i 

We first propose a dual picture for such observables. The Y-system will be constructed 
for the AdS^ case, with arbitrary number of operator insertions. In principle it should be 
straightforward to generalize to AdSc,, based on the result of the previous section. 

25 There is also an intriguing relation between the symbol of three-point form factor and six-gluon 
amplitude at two- loop at weak coupling 43 , 70] . It would be interesting to study this further at strong 
coupling, although this is not obvious by naively looking at the Y-system equations. 
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Figure 8: The dual momentum space configuration for the form factor with two-operator 
inserted. There are two periodic directions which give a periodic two dimensional lattice 
picture. This is equivalent to a torus. 

6.1 Evidence at weak coupling 

We first recall the picture of form factors with a single operator insertion. After T-duality, 
the picture involves a periodic null Wilson line boundary condition. The period is defined 
by the momentum of the operator. A duality between form factors and periodic Wilson 
lines was also found at weak coupling at one-loop [38J . A dual MHV rule description was 
proved for tree and one- loop form factors and also proposed to higher loops [ID] . 

How about a form factor with more than one operators? A natural generalization is 
that in the T-dual picture, every operator will generate a periodic direction. For operator 
Oi(qi), the corresponding period is qi. Such picture for the form factor with two-operator 
inserted is shown in Figure El It is given by a two dimensional periodic lattices and 
is associated to a Torus topology. For general m-operator insertions, the corresponding 
topology is T m . The momenta space can parametrized by introducing new coordinates 
x\ , which are defined as 

Xi% -Xi = q k , (6.2) 
l 

Q:=X>> (6.3) 
fc=i 

where k = 1, . . . , I, and I is the number of operators. See Figure El for the I = 2 case. 

One particular support of this dual picture is that the dual MHV rules description [iO] 
also applied to such generalized configuration, which provides an evidence that this dual 
picture may be applied more generally. 

In next subsection we will consider the worldsheet picture, where the introduction of 
multi-mono dor my seems to be a natural generalization of the single insertion case. As 
we will discuss later in subsection I6.5[ the monodromies in terms of the above spacetime 
coordinates can be also naturally defined. 



(fc (fc 
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6.2 Small solution and multi-monodromy 

We first introduce the picture of path 7^ on the worldsheet, as shown in Figure [9j The 
path 7^ is defined as a path going around the singular point z k where the operator Ok is 
inserted. A special path is which surrounds all poles. One may think this special path 
is similar to that of the single-operator case, ze^ 7o ° ~ ze l2lT , in the sense that effectively 
one can take the combination of all operators as one composite operator. 

Next we introduce small solutions Sj and s\ k . This is inspired by the spacetime picture 
of Xi and xf^ considered in last subsection, see Figure El 

The set of small solutions Sj is related to the special path 7^. As mentioned above, 
they behave similarly as those of the form factor with a single operator inserted, as one 
can take all operators effectively as a single operator. Therefore one has the same relations 

Sfl (z,() = B(()s (ze-^°°,O, (6-4) 

(^)M^,C) = n(C)(*)(*,C), (6-5) 

as the form factor studied in section |3j Using this set of small solutions, one can also 
define the T/Y-functions and construct the Y-system equations in exactly the same way. 
The monodromy Q should correspond to the product of the monodromies of all operators. 
The small solutions s[ is similar to the dual coordinate x\ in spacetime. Because 

(k) 

of the periodic structure, each set of small solutions s\ with fixed k is not different from 
the set of small solutions Sj. One has (s\ , Sj) = (sj, Sj). Similarly the Z 2 automorphism 
relations also applies such that sS(C) = i&3 s l ( e%n C)- The information of monodromies 
is encoded in the relation between Si and sj fe \ 

Recall that by definition, s^ + j is the small solution in the same sector as Sj but af- 
ter going around the complex z-plane (more exactly the path e^ 7 °° which surrounds all 
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operators) once, therefore they should be proportional with each other as given in ( 16.41) . 
Similarly, s^h is defined as a small solution in the same sector as but after going 
around the path 7^ once, as shown in Figure |9j One also has the proportional relation 
s h+i( z >0 s i{ ze ~ ^ lk i 0- O ne can introduce the proportional constant B^ as 

4 fc) (z,C) = BMsoize- (6.6) 

where k — 1,2, ... ,1. 

Because of the insertion of operators, small solutions are not single valued. For each 
path 7^, one can define a corresponding monodromy matrix as 

fyizefnt) = n<*>(C)(^)(*,C), (6-7) 

which is similar to the single insertion case. 
At the same worldsheet point, one has 

flf) (*,C) = B {k \0 ■ ^ k) )-\0 (z,C) , (6.8) 

B [k \0 = ( (5( y (C) B W (C) ). det[flW(C)] = l, (6.9) 

and similar relation between the s i+ n and Sj as (13. 9p . Since e^ 7o ° = nL=i e ^ 7fc > one h as 

1 1 
B = Y[ = JJ fiW. (6.10) 

k=l k=l 

Note that the order of operators should be not important, which implies that Q^ k ' should 
commute with each other. This requirement imposes further constraints on the mon- 
odromy matrices as will be discussed later. 



where 



6.3 New T and Y functions 

Now we consider the definition of T and Y functions and their relations. As we mentioned 

(k) 

before, if one just focuses on the set of small solution Sj and neglect s\ , one obtains a 
Y-system exactly the same as single operator case, with the total monodromy Q. The 
same Y-system can be constructed for the set of small solutions s\ with fixed k, since 
(si, Sj) = (s\ k \ Sj). The new degrees of freedom due to multiple operators insertions are 

contained in the interplay between small solutions Sj and s\ k \ which would necessarily 
involve the monodromies fi^. 
We define new T functions as 

rp(k) / (fc) \ rp(k) ; M\ + 

i l,2m+l :— \ S -m~l, S m+1 ), I l,2m : ~ \ S -m-l! Sm / , 

^~0,2m := ( s -m-l) s -m)i ^0,2m+l := ( s -m-2, + , (6-H) 

■— /« \ T (fc) ■— /« e (fc) \ + 

^2,2™, - — \*m)* m +l/) - t 2,2m+l - — \*m> *m+l/ ■ 
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Note that Tq^, T 2 \^, T[ are not normalized to be 1, and = (sq, Sq) ^ 0. Using Z 2 
automorphism, one also has the shifting relation (s i+ i, sS 2 ) = (sj, s^)^. 

Despite the difference, one still has similar Hirota equations by using Schouten identity 
(see appendix 1X1) 

rp(k)+ rp(k)- _ rp(k) rp{k) \TT ( Pi 1 9^ 

1 »,m -^ra - 1 o,m— 1 - 1 a,m+l - 1 a— l,m - 1 a+l,m • l u - xz 7 

Similarly, K functions can be defined as 

rp{k) rp{k) 

= T Q ' m -y m+1 ■ (6-13) 

a— l,m-L a+l,m 

The Hirota equations can be translated to the Y-system equations in a similar form as 

Y (k)+ Y {k)- /-, , y(fe) \/i ,y(fc) N 

'"./II ' -ll.l// _ 1 I / . i l ; — J 'd./// — .1/ (6 14) 



^o-l,m^a+l,m (1 + Ya-l,m)(l + ^a+l,m) 

In the normalization (sj, Sj+i) = (sj fc \ sS) = 1, the above equations are simplified as 

Yi k)+ Y^~ = (l + ^Oa + ^i), (6.15) 

where := y$>. 

6.4 Truncations and Y-system 

The main challenge is to construct a finite integrable system. We 'will first consider the 
truncation of Hirota equations, and then show how to write them into a gauge invariant 
Y-system. 

Similar to the single- insertion case (I3.12p . one can introduce the trace relation 

Tr[fiW(C)] = 5 (fe) (C) if } (C) - (B^)- l (0 Ti%(C) , (6.16) 

where ( = e -*( n + 1 ) 7r / 2 . can be expressed in terms of T^} 2 and Tr[f^ fc )], which provides 
a truncation for the chain of Hirota equations from the right-hand side. 
One can then define 

¥ {k \C) := (B^)-\C)Ti%(0, (6.17) 

and obtain 



7f(0 = (5 (&) )- 1 (C) [Tr[0«](C) + y W (C)J , (6.18) 

Y«\ = Tr[fi«]F (fc) + [F (fe) ] 2 , (6.19) 
where f2 (fc) (C) := ^ k \(). From this one has the equations 

y«+y«- = (l + Yi!))[l + Tr(^))F W + (F {fc) ) 2 ], (6.20) 

F ( fe)+F ( fc )- = 1+ y« 2 . (6 . 21) 
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Naively, one may introduce Y~® 2 and Y k for each new insertion of operators. However, 

the equation for YPj 2 contains Y^_ 3 , whose equation would then involve Y~<_ 4 and so on. 
This means that we also need to "truncate" the equations from the left-hand side, so that 
to formulate the equations into a finite system. 
We find it convenient to apply the relation 

/ S1 ,(*) ,(*) \ - /„ S W\ ( s ( fc > S W \ _ / Sl s ( fc ) w s M S W\ (a 22) 

\ s b s m-l/ \ s m > s m+l/ — \ a l!*m / \ a m-l! s m+l/ \ a l!*m+l/ \ 3 m-l) s m / ) 

which written in terms of T-functions is (up to a phase shift) 

rp(k) _ rp{k)+ rp[m+3] _ rp(k)[2] (a r,q\ 

1 m — 1 m+l 1 l 1 m+2 ■ {V.Z,0) 

This provides a recursion relation for Tm \ and one can express all in terms of only 
two functions and 7\. For our purpose it is enough to consider!^! 

4% = Tf 1 - rf_f . (6.24) 

Together with the trace relation involving one can truncate the chain of Hirota 

ffc) (k) 
equations involving only T^_ 2 and T~_ v 

We need to further write the truncated Hirota system into a gauge invariant Y-system. 
To do this, one can first write Y^_ 2 as 

v (k) _ T (fe) T {k) _ (T (k)+ T [h\ _ T (fc)[2]x T (fc) 

1 h-2 ~ I n-3 I n-l ~~ K 1 h-2 ± l 1 h-1 ) 1 n-1 

= Z ^ - 1 - Y^+ , (6.25) 

where a new function Z^ is introduced as 

Z ( k ) T W+ T W iT f. ( 6 .26) 

The advantage of introducing Z^ is that it is straightforward to obtain the equation 

Z (k)+ Z (k)- = ^(k) _ yikj + y (y(k)^ + ^ ^[n] + ^ {6 27) 

In particular, yPj 3 no longer appears, and one obtain a closed set of equations. Therefore, 
rather than use yPj 2 , we will use the new function Z^ k \ 



Y {k) - 



To summarize, one obtains a closed Y-system with the functions Y m , Y, and Z^ k \ 



Y+Y- = (l + Y m+1 )(l + Y m _x), m = l,..,n-3, (6.28) 
^- 2 *7- 2 = (l + Yi-3)(l + Tr[fi]Y + Y 2 ), (6.29) 
Y + Y~ = 1+Y ft _ 2 , (6.30) 



9fi (k) (k) (k) 

b One may also use this relation to solve for Tl in terms of Tl_ 2 and and then substitute it 

into Y^_ x = Th-2Th- This will give the same equation as (|6 . 25[) . 
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and 



Z (k)+ Z (k)- = ^ Z (k)+ _ yW[2]^ (y(k^ + ^ (y[n] + ^ ^ ( g 

- (fe ) + -( fc) - = z (fc)_ y W+ ? ( 632 ) 



where 



= Tr[fi fc ]r (fe) + [r (fc) ] 2 , (6.33) 

and k = 1, . . . ,1 — 1. One can see the function in (16.3ip has phase shift which is 
beyond the physical strip (— ir/2, 7r/2). Therefore to write it into a integral equation, the 
extra pole contribution should be considered. 

We comment on the degree of freedom. One can see that for each operator Ok, two 
new functions are introduced. One may understood this using the same argument of the 
single insertion case: the new 2x2 unitary matrix subtracting the trace Tr [Q^] leaves 
two independent components. However, there are also extra constraints: the monodromy 
matrices should commute with each other. This implies the matrices should be in general 
in the form of 

V o 

b k l/a k/ 

and only one new degree of freedom is introduced for each new operator (for example, 
given bk then a k is fixed). This implies that the Y k and are actually not independent. 

This matches with the degrees of freedom from spacetime boundary configuration 
that we considered in subsection I6.1[ where each new operator introduces a new periodic 
direction characterized by q. Note that one F-function in AdS$ gives to two real degrees 
of freedom, due to the left and right hand decomposition. Since the boundary information 
enters into the Y-system via the WKB approximation, this is also related to the structure 
of P(z) functions which will be discussed in section [7J 



(6.34) 



6.5 Spacetime picture 

Now we consider the monodromies in terms of spacetime variables. We first recall the 
dual momenta space configuration 

_ _ (*0 _ ( fc ) _ (*0 _ _ (a 9c\ 

%i — x i — Pi ) x i+n X i — Qk > ^O.OOJ 

I 

Q = $>. (6.36) 
i=i 

Each monodromy corresponds to a conformal transformation which maps {X^\ X^li, X^ 2 } 
to {Xq, Xi, X2} and can be defined in terms of left-hand spinors as 

\(k),L A(fc) X L . ( &11>%1 + ^U^iM / fi o 7 \ 

Vhi \a + "22 \,2 



(k) (k) 
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where i = 0,1,2. One has 

*•-&!)■ fl -6 5- »*> 

which are indeed commuted with each other and satisfy Q = nf=i 

One can express and in terms of spacetime variables which specify the shape 
of dual Wilson line configuration. For the Y k functions, it is similar to Y 

/i \ l\L + _ (fc) + 

y«( C = i A+i) = - } Al ' A \ l = - I' % m+ . (6.39) 

For the function Z^ fc \ one can first write it in the gauge invariant form 

T (k)+ T (k) T [h] , q ( k h/ q »( fc )\/,( fc ) \ 
zW(c) = = ( a o,* fl )( !' A J ( ^ r (Ce-' 4 ). (6-40) 

rp rp rp[n\ I \ / ( ) ( fc )\ / \ 

J0,n-2 ^2,n-2 ^ 2,0 \ S 0, Si/ S ri / \ S n ' S n+l/ 

Then the spacetime expression can be given as 

(x - Xi) (4^! - 4 7 ) (4 - X^J 
which is manifestly conformal invariant. 



7 Function P(z) and WKB approximation 

As reviewed in section the boundary conditions are related to the holomorphic function 
P(z) which is also related to the WKB approximation. In this section, we study this in 
more details. We will focus on the cases with short operators, which are dual to periodic 
Wilson line configurations. The general structure of P(z) will be given. The general 
pattern of corresponding WKB lines will also be studied. 

7.1 P(z) for general form factors 

For amplitudes or null Wilson loops, P(z) is a pure polynomial. For form factors, due 
to the insertion of operators, new pole terms are involved. This can be understood by 
studying the behavior of the solution near the horizon. Our discussion is for general AdS^, 
following the AdS^ case in [36J. 

The main trick is that by doing a worldsheet conformal transformation, one can bring 
the horizon at infinity to the origin in the new coordinate. Firstly we recall the picture 
of the dual surface in Figure [Q Without loss of generality, one can take the periodic 
direction to be along x\. We parametrize the worldsheet by the coordinate z = r + ix\. 
Near the horizon where r — > oo or z — > oo, the surface is asymptotically a straight strip. 
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We can set x-i = £3 = t = up to a translation. The induced Poincare metric takes the 
form 

dsL = j0^- 2 • (7.1) 

For this simple solution, the corresponding polynomial is simply P(z) = 0. 

One can now apply a standard coordinate transformation to map the strip to the unit 
disc with new coordinate z 

z = e~~ z . (7.2) 

The infinite of z becomes the origin of z. It is in this new coordinate z that we discuss 
the picture of the worldsheet monodromies in previous sectionsEll- In the new coordinate, 
the above induced metic takes the form 

rfs i 2 nd = :r~T7 — \ = e~ 2a dzdz , a = - ^ \og(zz log 2 (zz)) . (7.3) 
zzlog (zz) I 

This provides the boundary condition for the solution when z — > 0. In particular, one can 
see that unlike amplitudes, a is not regular near the origin, which means the worldsheet 
is no longer smooth. Near the cusps when z — > 00, one still has a — > 0. 

In AdSs, the function p(z) ~ iV ■ d 2 X |24J and have the transformation property 

V{z) = (t^ p(z) = z 2 p{z). (7.4) 



dz J 

As discussed above when z — > 00 (or z — > 0), p(z) — > 0, this implies that 

p( z ) = - + O(z ). (7.5) 

z 

For the AdS's case, P(z) = d 2 X ■ d 2 X which gives 



dz 



-4 



P(z) = ^-j P(z) = z*P(z). (7.6) 

The condition P(z) = when z — > 00 requires that 

P(«) = ^ + ^ + ^ + 0(z°). (7.7) 

However, 1/z 3 term is not allowed. This can be understood as that near horizon, the 
AdS 5 solution can be embedded into AdS 3 , which must then satisfy P(z) oc p{z) 2 [25] . 
One concludes that for form factors in AdS^, P{z) has the following general structure 

P(z) = a n ^z n ~ 4 + • • • + a x z + a + — + \ . (7.8) 

z z z 



27 For small solutions the boundary is at \z\ — > 00. It seems there would be a problem as \z\ — > 00 
implys \z\ = r — > — 00. As the focus here is on the behavior near the horizon, one should think that above 
transformation is only for the region near the horizon. 
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For n = 2 case (which can be always embedded in AdS^), one has 




P {zf 



1 



(7.9) 



z 



For the three-point case it is given as 




a_i 



1 



(7.10) 



+ 



z 



z- 



For a general n-point form factor, the number of free parameters from the coefficients 
cii is 



For the AdS± case this matches exactly with the degrees of freedom from a counting of the 
symmetries of a periodic null Wilson line configuration. For AdS§, there should be further 
(n — 3) parameters from gauge connection, as in the case of scattering amplitudes • In 
total the number of the parameters is 3n — 7, which is also consistent with the counting 
of symmetries. 

For multi-operator insertions, a natural proposal is that each operator introduces one 
new pole term. For example, for the case with two cusps and two operators we would 
have 



* z ~ z o 

where we have used scaling and translation symmetries on the worldsheet to set = 1 
and Zq 1 ^ = 0. In the limit that z^ — > 0, it reduces to the single operator form, which is 
consistent with the picture in section |6j 

In this proposal, the degree of the polynomial is related to the number of cusps, and 
the number of poles corresponds to the number of operators inserted. Each P(z) function 
define an algebraic curve, or a Riemann surface. The numbers of genera and singularities 
are related to the numbers of cusps and operators. It also produces a consistent WKB 
line picture as shown in the example in next subsection. 

However, this seems to be not the full story. The problem is that the remaining two 
complex parameters in (I7.12p give four degrees of freedom, which do not match with the 
T-dual spacetime picture which has only 2 degrees of freedom. This implies that one may 
need to impose extra constraints on the coefficients a_\ and Zq for each insertion. This 
seems to require a better knowledge of the T-dual picture of the minimal surface from 
which one may do a similar study as for the single insertion case. 

There is also another possibility. Although the function P(z) contains more parame- 
ters, the final area may be independent of these extra degrees of freedom. In other words, 
some of the parameters in P(z) may be taken as "gauge-like" degrees of freedom, and one 
may change them without changing the physical area. This picture seems more natural 
but needs to be checked through a detailed study of the area. 

In either case, we believe that the general structure of the P(z) function should be 
correct. We summarize the function P(z) for general form factors in Table [TJ 



2(n - 3) + 2 . 



(7.11) 
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AdS 3 (n = n/2) 


Amplitudes 


p(z) = Z h ~ 2 + CLn-4, 2™" 4 H h Oi z + a 


One-operator 


p( z ) = z h ~ 2 + a h _ A z-- 4 + • • ■ + a + ±± 


Multi-operator 


p( z ) = z *-z + a^_ 4 z"~ 4 + ■ ■ • + a + E,-i -^ll 




AdS$ and AdS± 


Amplitudes 


P( z ) = z n ~ 4 + a n _ 6 z n ~ 6 + --- + a lZ + a 


One-operator 


P(z) = z- 4 + a„_ 6 z^ + ■ ■ • + a + M + ^ 


Multi-operator 


P(z) = + a n _ 6 + ■ ■ • + a + Er=i 



Table 1: JTie P(^) function for amplitudes and form factors in AdS 3 and AdS 5 . n is the 
number of cusps, m is the number of operators inserted. 



7.2 WKB approximation 

The asymptotic behavior of Y-functions is determined by P(z) through WKB approxi- 
mation. The WKB lines are defined by the parametric line z(t) as: 

where 

AdS 3 : p{z) = x 2 , AdS 5 : P(z) = x 4 . (7.14) 

The AdS 3 case corresponds to the SL(2) Hitchin system which was studied in details 
in [57]. As 9 changes (( = e 10 ), the WKB lines will change correspondingly which is 
related to the wall-crossing phenomenon in M = 2 theory. Although the physical context 
looks quite different here, the mathematics is basically the same. Below we summarize 
the general patterns of WKB lines for both AdS 3 and AdS$ cases. 

The WKB lines in Ad S3 have the following structure. For a general point, there is 
only one WKB line going through. The special points are the zeros and poles. There 
are three lines ending on each zero, and one line ending on each (simple) pole. These are 
shown in Figure [10J The full WKB lines for the six-point form factor with two-operator 
inserted is shown in Figure [TTJ 

The AdS§ case corresponds to a SU(A) Hitchin system. The WKB lines have more 
complicated structures. For a general point, there are two WKB lines going through^ 



3 To be more precise, this is the picture projected on a single z-plane. P(z) — x defines a Riemann 
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Figure 10: WKB lines in the AdS^ case, where we choose 9 = tt/2. The two figures 
illustrate the behavior of WKB lines near zero and simple pole respectively. WKB lines 
which end on zeros or poles are shown in orange color. Zeros are shown as blue points, 
while poles are denoted as red points. 




Figure 11: WKB lines in AdS^ for the polynomial p{z)= ^fez^r^ — -, = ir/2. It 
corresponds to a six-point (n=3) form factor with two operators inserted. The second 
figure plots the WKB lines in a much larger range, which makes it obvious that there are 
three small solutions at infinity. 



There are five lines ending on each zero, three lines ending on each simple pole and two 
lines ending on each double pole. The WKB patterns are show in Figure O The WKB 
lines for a four-point form factor are given in Figure [T3J 

One can associate small solutions to to the asymptotic WKB lines, as labeled in Figure 
[TT1 and [121 Integrals over WKB lines (show as black lines in the figures) that connect 
different small solutions will provide the WKB approximation for the contraction of the 
small solutions. For Y-functions, the corresponding WKB lines always form a closed 
contour. Therefore, the WKB approximations of Y functions in the limit of ( — > 0, oo are 
given as cycle integrals, which is related to the mass parameters [26] . They are related to 
the coefficients a, appearing in P(z), and also implicitly related to the shape and periods 
of Wilson lines. 

surface with four branch covers. On each sheet there is only one WKB line going through each point. Four 
sheets give actually four lines. Two of them overlap with the other two (but with different orientations) 
[26) . Projectively one gets the figures shown here. Similar picture applies for the AdS^ case. 
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Figure 12: WKB lines in the AdS$ case, where we choose 6 = 0. The three figures 
illustrate the behavior of WKB lines near zero, simple pole and double pole respectively. 
Orange lines are special WKB lines which end on zeros or poles. 




Figure 13: WKB lines in AdS$ for the polynomial p(z)= ^ z ~ 1 'i z ~ 2 ' , 6 = 0. It corresponds 
to a four-point form factor. The second figure plots the WKB lines in a much larger 
range, which makes it obvious that there are four small solutions at infinity. 

8 Summary and discussions 

In this paper, we study form factors in M = 4 SYM at strong coupling in AdS& and with 
multi-operator insertions. These are two non-trivial generalizations of the AdS% form 
factors studied in [36] . 

The generalization to AdS$ involves new technical problems comparing the simple 
AdS% case. The main challenge is how to introduce the truncation conditions with a more 
complicated 4x4 monodromy matrix, and how to write the system in a gauge invariant 
way, i.e. in terms of Y-functions. We clarify and solve these problems. The Y-system of 
three-point case is constructed explicitly, which potentially would have a connection to 
the strong coupling Higgs-to-3-gluon amplitudes in QCD. 

The second generalization to the multi-operator insertion cases provides a more inter- 
esting picture and we would like to make a few further comments. The main hope is to 
provide a new technique to study correlation functions. In doing so, we take an uncon- 
ventional point of view at strong coupling: to apply the on-shell techniques to compute 
off-shell observables. Although similar ideas have been studied in the weak coupling side, 
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this point of view seems have not been taken seriously at strong coupling. 

According to the picture we proposed, adding one operator corresponds to introducing 
a new monodromy matrix, which is taken as a condition imposed on small solutions defined 
by on-shell cusps. The techniques developed for amplitudes or null Wilson loops can be 
applied to computing such more general class of observables. For AdS^, we construct 
the Y-system explicitly for the cases with arbitrary number of operator-insertions. The 
construction should in principle be generalizable to AdS§ based on the prescription of the 
single-operator result developed in this paper. 

The derivation of the Y-system should be in principle applicable for general operators. 
Different operators would correspond to different kinds of monodromies. The simplest 
case is for short operators which are dual to light string states as having been studied. 
For them, the monodromy depends only on the momenta of the operators. This is actu- 
ally very interesting, considering that normally it is hard to study correlation functions 
with purely light operators at strong coupling besides the perturbative Witten diagram 
techniques, partially due the complexity of string vertex operators in AdS$ (see some 
ideas in |71j). It would also be interesting to construct the monodromy for more general 
operators, in particular classical solutions such as the GKP string [72]. 

Although our construction relies on the partially on-shell structure of the observables, 
the multi-operator structure should in principle contain all kinds of information in corre- 
lation functions. In particular from the general structure of OPE 

1 (x)0 2 (y) = ^c 12a f(x-y,d y )O a (y), ( 8 - X ) 

a 

an immediate step would be extracting the OPE coefficients by using form factors with 
two operator insertions and comparing it with form factor with the single operator O a in 
the OPE expansion. 

Although our construction of Y system does not rely on an exact knowledge of the 
string solutions, we do not have an explicit string T-dual picture of a form factor with 
multi-operator insertions. As being discussed in section [7J this should provide a better 
knowledge for the function P(z). The T-dual picture of short string states was studied 
recently in [73], where a similar Wilson line picture was obtained. It would be interesting 
to understand the picture of interacting multi-closed-string states. 

While there are lots of studies on correlation functions, we would like to point out 
particularly [7jJ[7j2[76j[77], wnere quite similar integrability techniques have been used. 
However, the detailed physical pictures and the building blocks are quite different. The 
method in those papers is limited to classical heavy operators (where the S 5 geometry 
plays also an important role), while our prescription is focused on short operators (but 
in principle could be for more general operators). It would be interesting to study the 
possible connection between the two pictures. 

There are interesting algebraic curves appearing in the construction as discussed in 
Section [7J In gauge theories similar algebraic curves (the Seiberg-Witten curves) also 
appear [67]. It would be interesting to study their possible connections. There are also 
similar spectral curves for classical string solutions such as those studied in [78]. There 
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is an important difference though: while the spectral curve is a curve defined on the 
spectral parameter £-plane, the algebraic curve here is on the worldsheet z-plane. On the 
other hand, in our picture there are close interplays between the two planes. It would be 
interesting to understand the connection more explicitly. 

We would also like to make a few comments on the symmetries. Unlike amplitudes, 
form factors do not have dual conformal symmetry^. However, as we have seen that there 
is no problem to use integrability to compute strong coupling form factors. Technically this 
may be understood that through changing of variables, one can bring spacetime quantities 
to a picture on the worldsheet, where the symmetries are in some sense enhanced and 
integrability techniques can be applied. It would be very interesting to study its possible 
correspondence at weak coupling, for example to have a realization of small solution 
picture at weak coupling, see an interesting proposal along this direction in In our 

opinion, it would be the symmetry of the theory rather than the symmetry of observables 
that plays the most important role. 

Finally, let us mention that there are a few technical problems to clarify. While we have 
obtained the Y-system and considered the WKB approximation, we have not discussed 
how to find the explicit solutions. The explicit integral form of the Y-system equations 
are not given. The main complexity is due to the phases of some functions appearing 
in the equation are outside the physical strip, and extra pole contributions need to be 
considered. It is also necessary to show how to derive the expression for the area. For the 
AdSs form factors, a natural expression for the non-trivial free energy part is proposed, 
while for the multi-operator case a further study is necessary. There are also some issue 
about the P(z) function of the multi-operator insertion case. These problems are under 
investigation and we hope to report them in the near future. 
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A Hirota and Y-system equations 

In this appendix we review the definition of T and Y function, and some basic equations 
among them. 

29 This may be understood most easily in the T-dual picture, where a periodic Wilson line does not 
preserve special conformal symmetry. 

30 Scc also [80 for an interesting idea of introducing spectral parameters for amplitudes at weak coupling. 
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We start with Cramers rule: 

S h ( S «2> S «3> > S *fc+l) _ S «2 5 S «3> 5 S ik+l) + ' ' ' + ( — 1) S i k +1 ( S il> 5 S «fe) = ) (A.l) 

where s, is a dimensional vector and the contraction is defined as 

(s 1 ,s 2 ,. . . ,s k ) := e Ql -° fc si iai . . . s k<ak . (A.2) 

Pliicker relations can be obtained by contracting the small solutions with another set of 
small solutions s^, ■ ■ ■ , Sj k _ 1 

(Sj 1 , ■ ■ • , Sj k _ 1 , Si 1 )(Si 2 , Si 3 , • ■ ■ , Si k+1 ) + ■ ■ • + ( — 1) (Sjj, ■ • • , Sj k _ 1 , Si k+1 )(Si 1 , ■ ■ • , Si k ) = . 

(A.3) 

When k = 2, one gets the Schouten identity 

(si,Sj)(s k ,si) + (si,s k )(si,Sj) + (si,si){sj,s k ) = 0. (A.4) 
When k = 4, one obtains useful relations for the AdS§ case, such as the Wronskian relation 

\Sj, Sj , S a , Sft) (Sfc, S; , S a , Sfo) -|- (Sj , Sfc, S a , Sft) (s; , , S a , Sfc) -|- (Sj , Si , S a , Sft) (sj , Sfc, S a , Sfr) . 

(A.5) 

Below we give the definition of T- and Y-functions, then we apply the above relations to 
obtain corresponding equations. 

A.l The AdSz case 

We use the convention (note it is different from the AdS§ case): 

/± := /(e^C), f [k] := f(e^C) • (A.6) 
In this notation one has from (I2.19P 

(s i+1 ,s j+1 ) = (s i ,s j } [2] . (A. 7) 

The T-functions are defined as 

7l,2m+l := ( s -m-l; s m+l)) Ti 2m := ( s -m-l> s m) + , 

To t 2m := ( s -m-l> s -m)i ^0,2m+l := (S-m-2; + ) (A. 8) 

^2, 2m := ( s m 5 s m+l)) ^2,2m+l := ( s m) s m+l) + - 

Ti jm is non-zero for m = 1, ...,n — 1, and the normalization (sq, si) = 1 corresponds to a 
gauge choice T 0jm = T 2 , m = 1. 

Using Schouten identity (1A.4[) . one can obtain the so-called Hirota equations 

^ a,rr\F a,m = ^a,m— l^o,m+l + ^o— l,m^a+l,m > (A -9) 

where the indices a, m take integer values. 
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Hirota equations contain huge gauge redundancies 

TaAQ -> U 9^ {aa+ " m) C)T a , m (0: (A.10) 

a,/3=± 

where g a p(C) are f° ur arbitrary functions. Like defining field strength in gauge theory, 
one can introduce gauge invariant functions, so-called Y-functions 

T am — 



Y 

± a,m rri rri 

■L a—l,m-L a+l,m 

The Hirota equations become the equations for Y-functions 

,m— 1 ,m+l J 



y+ y- 



(A.11) 



(A.12) 



Y a — l,m^o+l,m (1 "I - ^a— l,m)(l ~1" ^o+l,m) 

In the normalization (sj,s i+ i) = 1, Yo, m ,Y2, m are trivial (either zero or infinity), the 
equations of Y-functions simplify as 



Y+Y~ 

m m 



(i + y m _x)(i + y r 



m+lj j 



where Y n 



(A.13) 



A. 2 The AdS 5 case 

We use the convention: 

/ ± := /(e^C), / [fc] := /(e^C) ■ 
Useful relations due to the Z 4 automorphism are 

(Sj-i, Sj, Sk-i, Sfc)^ = (Sj, Sj+i, Sk, Sfc+i) 
(s fc ,S J _2,Sj-l,Sj) [21 = (Sfc, Sfc+i, Sfc+2, Sj) 

Define the T functions as 



Ti, m (C) 

^2,m(C) 
T 3 ,m(0 
Ti,m(t) 

Using Pliicker relations, one obtains the Hirota equations 
Gauge invariant Y-functions can be defined similarly as 



(8_ 2)5 _ 1)5o ,Sl) [ - m - 11 ) 
(S_ 2 , S_i, S , S m+ i) [ ~ m] , 
(S-l, So, S m+ i, S m +2)' ™ 1 ' i 
(S— 1) S m , S m _j_i, Sm+2)^ ^ 5 
(Sm; S m -)-l, S m _|_2, Sm+3)^ ^ 



a'm 4—a,m 



1,2,3. 



y ■= 



T T 

ajm+l- 1 - 4— a. 



m— 1 



The Hirota equations become the Y-system equations: 

a + Y a 

,m+l )(l + ^4 —a,m—l ) 



y- y+ 

a,m A:— a'm 



Y a +l,mY a - 



l,m 



'1 + yo+i,m)(i + y-i,m) 



a = 1,2,3 



(A.14) 

(A.15) 
(A.16) 



(A.17) 



(A.18) 
(A.19) 

(A.20) 
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B Twistor variables 



In this appendix we given a brief review on (momentum) twistor variables, see for example 
[62JI8T] . One technical point we would like to clarify is how to transform twistor variables 
into Lorentz variables, and vice verse. 

We first recall the relations between embedding and Poincare coordinates 

X» = — , X + = -, X- = , (B.l) 

where (77^ = (-1, 1, 1, 1)) 

- 1 = -X+X- + X"X„ , X ± := X- 1 ± X 4 . (B.2) 

Twistor variables can be understood as in the spinor representation of the embedding 
5*0(2,4) space, which is a S77(4) fundamental representation, denoted by A, 

With one explicit choice of gamma matrices, one has 

/O iX + X 2 + iX 3 X 1 -X° \ 

1 

* * iX~ 

\* * * J 

Note that det(X) = (X.X) 2 /4. Due to the freedom of choosing normalization, twistor 
variables Aj are projective coordinates in CP 3 . 

Using (IB. II) and (1B.4I) . one can obtain the relations between A a and x^. For example, 
define 

x ± = x 1 ± x ° , (B.5) 

one has the relations 



Xab ~ 71 



A [0 A 6] . (B.4) 



. X iM . A[ X A^ 1 Xi 23 . A [2 ^3]~ 1 m „ N 

X- = 2— =1— r— Xj = I— =%—. r— . (B.O) 

Another description of so-called momentum twistors is practically more useful which 
was first introduced at weak coupling [62] Ell- Consider a null Wilson line configuration 
defined in the momentum space of amplitudes or form factors 

x i+l ~~ %i = Pi 5 Pi = j Pi,a6t = Aj j0 A^q, , (B-7) 

where the left and right-hand SU(2) Weyl spinors are denoted by A and A. We define 
Xij := Xi — Xj. The Weyl spinor contractions are defined as 

:= e Q/3 A ija A ii/3 , := e^A^A-^, (B.8) 

(i\xy\j) := A? ./v,//".^.,-. := A> a/ jA*. (B.9) 



31 It is called momentum twistor just because it is denned in the momentum space, mathematically 
there is no difference to usual twistor. 
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The momentum twistors can be explicitly denned as follows 

\i = (A iia , fj, i)& ), Hi,a ■= -i(xi-Ai)a = -ie af3 Xi t0l aA it p . (B.10) 
Note that also /i^ = —% (xj+i • Aj)^ . The contraction of twistors is defined as 

(Aj, \j, Afc, Xi) := e abcd \ i: a\j tb \ k:C \i !d , a=(a,a). (B-ll) 

The geometric picture of twistor space is that, each spacetime point corresponds to 
a line in twistor space determined by two twistor variable, X{ ~ Xi ~ Aj_i A Aj. If two 
spacetime points are null separated, the corresponding two lines in twistor space intersect 
with each other. This is obvious in the above definition since null-separated Xi and Xi + i 
both contain Aj. 

To write the contractions of twistor variables in terms of Lorentz coordinate, one 
practically very useful formula is [81] 



(i\x i>jXj>k \k) = {Xi > y- 1 ;^; Afc) . (B.i2) 

For example, using flB.12j) it is easy to obtain the relation 

- = (K-iAiAj-uXj) (B13) 



Furthermore, any normalization independent expression of twistor contractions can be 
written in terms of Lorentz variables. For example, for the ratio variables appearing in 
form factors, one has 



(Ai, A 2 , A 3 ,^A 4 ) (Ai, A 2 , A 3 , A 4+n ) (l|p 2 (g+_P3)j4) 

(Ai, A 2 , A 3 , A 4 ) 64 (Ai, A 2 , A 3 , A 4 ) (l|p 2 P3|4) 

(Q Ai, A 2 , A 3 , A 4 ) (A n+i , A 2 , A 3 , A 4 ) + p i2 ) p 3 |4) 

(Ai, A 2 , A 3 , A 4 ) h (Ai, A 2 , A 3 , A 4 ) (l|pi 2 p 3 |4) 



(B.14) 
(B.15) 



where we use the relation (I4.55P X i+n = biQ Xi. Using flB.12j) it is easy to perform the 
computation for the middle expressions. 



C Monodromy with a different basis 



In this section we briefly explain the monodromy defined in a different basis of small 
solutions, in particular how Q is related to fl. 
First recall the definition of the monodromy 

/ s i \ / s i \ / s 2 \ / s 2 \ 

So (ze 2m ,() = n- 1 (C) ( :. c. ) . 

s -i 

\s-2J 



(8i\ 




(s 2 \ 


so 




Si 


S-l 


so 


\s-2J 




\s-i) 







(H [2] )- 



\0 



Sl 

so 
\s-xj 



(C.l) 
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and 



(s n+ i, s n , s n _i, s n _ 2 ) («,C) = B X (C) (si,s ,s_i,s_ 2 ) (-2e 27rl ,C), 
(sn+2,s n+ i,s n ,s n „ 1 ) T (z,C) = (B 121 )" 1 ^) (s 2 ,si,s ,s_i) T (ze- 27ri ,C) , (C.2) 



where using the relations of (I2.20p . the proportional constants are given by a single B 
function 



B- 1 = &&g{B,{BW)-\B^\{B^)- 1 } 
B~ l = tiag{(BW)-\Bl- 2 \B-\B®}. 

One can expand s_ 2 in terms of s , s\, s 2 }, then one gets 

S-2 — T X1 s_i — T 2j i Sq + t\ x si — s 2 . 



(C.3) 



(C.4) 



Similarly for s n _ 2 one has the expansion 

) Sn-1 ~~ ^2,1 S ™ + (Sn-2, S n -1, S n , S n+2 ) S n +i — S n+2 , (C.5) 

where the other two contractions are Ti i or T^i depending on whether n is even or odd. 
By introducing 
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